Kob-Andersen model: a non-standard mechanism for the glassy transition 
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We present new results reflecting the analogies between the Kob-Andersen model and other glassy 
systems. Studying the stability of the blocked configurations above and below the transition we also 
give arguments that supports their relevance for the glassy behaviour of the model. However we 
find, surprisingly, that the organization of the phase space of the system is different from the well 
known organization of other mean-field spin glasses and structural glasses. 



I. INTRODUCTION 



In the last two decades our understanding of the glass 
transition in spin systems has provided a coherent picture 
of glassy phenomena at the mean field level §]. These 
systems provide the simplest example of systems where 
complex energy landscapes, with many local energy sta- 
tionary points, lead to vitrification and to the absence 
of thermalization over large time scales. At a mean field 
level, there are spin glass models which show a dynamic 
transition very well differentiated from the static one. 
Below the dynamical transition temperature the system 
is unable to thermalize and get stacked at energies above 
the equilibrium one. At the same value of the tempera- 
ture the space of equilibrium states gets disconnected. 

The situation is different for real materials, where the 
sharp transition is not present, and seems to be substi- 
tuted by a smooth crossover to a dynamics dominated 
by activated processes 0. Progresses in the last decade 
of research, however, indicate that many of the aspects 
of the picture of the glass transition gained from the 
study of mean field models may well apply to real ma- 
terials. In fact, while already Kirkpatrick et al. ||]4|] 
noticed the similarity between the thermodynamic tran- 
sition in structural glasses and spin glasses, later there 
were found discontinuous transitions in spin glasses with- 
out quenched disorder giving confidence that this anal- 
ogy was not fortuitous ■ Now, it has been recognized 
that at least in some approximation, classical first prin- 
ciple models have a phenomenology which resembles to 
that already found in mean field spin glasses |Q,0j8| . 

Moreover, an important feature of the phenomenology 
of glass forming liquids is the formation of dynamical het- 
erogeneities, characterized by a spatial extension and a 
given life-time, both of which increase as the temperature 
is decreased. These heterogeneities have been found in 
numerical simulations of various realistic models H[n| as 
well as in experiments on colloidal glasses |llj , and seem 
to be a salient feature of landscape dominated glassy sys- 



tems. In order to quantify the presence of dynamic het- 
erogeneity, it has been proposed to study 4-point density 
dynamic correlation function. For this quantity mean 
field theory predicts a peculiar time/temperature depen- 
dence, in agreement with the one observed in numerical 
simulations Jl^| . 

A related problem is the physics of granular systems. 
Since, in general, these systems involve many particles it 
is tempting to treat them using the methods of statistical 
mechanics and, in fact, this approach has become an ac- 
tive field of research for people with a background in this 
field [fl3|| . Particularly interesting have been the experi- 
mental results of the group of Chicago who showed that a 
granular system subject to continuous tapping develops 
a glass-like transition [fblf . 

An alternative view of the glass transition comes from 
kinetic models, which have been useful to prove that com- 
plex energy landscapes are not a necessary ingredient to 
have glassy phenomena. Vitrification has been found in 
simple lattice cellular automata models with kinetic con- 
straints. 

An interesting model belonging to this class is the Kob- 
Anderson model JT5| , which has been shown by several 
studies to undergo a structural arrest at densities larger 
than the threshold one. The phenomenology of the model 
has been compared with qualitative agreement with the 
Mode Coupling Theory which proposes an intrinsi- 
cally kinetic view of the glass transition. 

Despite the fact that the mechanisms leading to struc- 
tural arrest in the model can not be related to a non- 
trivial energy landscape, many aspects of the physics of 
the model are very similar to the phenomenology of mean 
field glassy systems. In particular the blockage to the 



threshold density found in aging experiments in [|16 17 
closely reminds the blockage of the internal energy in 
the dynamics of mean field spin glasses. In addition, it 
has been recently shown that during aging dynamics, the 
system probes portions of phase space close to generic 
random blocked configurations [[L8| , again, in a way rem- 
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iniscent of mean field glassy phenomenology. 

It is tempting in such situation to invoke some univer- 
sal mechanism underlying structural arrest both in land- 
scape dominated systems, and in kinetically constrained 
ones. In this paper we investigate the limits to which 
these analogies can be pushed and to what extent the 
mechanisms are really universal. 

The Kob-Andersen model, that first emerged as a 
model to describe the dynamic transition in structural 
glasses, has also been recently used as a simple model for 
the jamming transition in granular systems ]l9|Jl^^8[ . 
While many other lattice models have been proposed to 
study similar problems |2^,^l|, this one has the advan- 
tage of a trivial equilibrium state that can be easily char- 
acterized analytically without further assumptions, and 
that therefore allows to a clear comparison between the 
static and the dynamic properties of the system. 

A first test to which we submit the KA model is the 
computation of the dynamical susceptibility associated to 
a 4-point function, and show that it behaves as expected 
from mean field theory. 

A second point that we investigate is the nature of the 
blocked states. A lot of emphasis has been recently made 
about the role of stationary points of the energy surface 
in glassy dynamics, and the ideal mode coupling transi- 
tion has been associated to the point where the number of 
negative directions of the inherent saddles vanishes |^2|,0 ■ 

In other words, while above, the MC point the dynam- 
ics spends a lot of time close to saddles, below, it spends 
its time close to minima. The work of |Q has empha- 
sized the role of the blocked states during aging in the KA 
model. In this paper we show that while this states are 
responsible of the glass-like properties of the model, their 
geometrical organization differs from the organization of 
the states found in spin and structural glasses models. 

The remaining of the paper is organized as follows: In 
section [n] we define the model and discuss the current 



understanding in the literature. In section III we present 



our numerical simulations and discuss the results, finally 



in [V the conclusions arc outlined. 



II. THE MODEL 

The Kob-Andersen model is a lattice gas defined in a 
three dimensional lattice of size L 3 where N particles in- 
teract by hard core repulsion (i.e only one particle can 
be in a site at a given moment) and that evolves using 
dynamic rules (described below) which are symmetric in 
time, so that the detailed balance is satisfied. 

We first choose a particle at random, then the particle 
moves if the following conditions are satisfied 

• the neighboring site is empty 

• the particle has less than 4 nearest neighbors 



• the particle will have less than 4 nearest neighbors 
after it has moved 

Already in the first paper devoted to the model, it was 
shown that the dynamics becomes slower when the den- 
sity p in the system increases. In particular the diffusion 
coefficient vanishes as a power low at the critical density 
p c m 0.88 independently of the system size ruling out the 
possibility that this phenomena could be associated with 
the formation of an infinite backbone in the system [fl5| . 

More recently, the system has been studied in the pres- 
ence of a chemical potential [Q. Despite of the trivial 
equilibrium properties of the system, it was proved in 
a series of papers that: Below a certain critical density 
[p c w 0.88) the system is not ergodic, independently of 
the chemical potential the density gets blocked at a den- 
sity lower that the critical density. This density depends 
on the "cooling rate". Detailed linear response numer- 
ical experiment have shown the emergence of an effec- 
tive temperature. This effective temperature coincides 
very well with the "Edwards" temperature, which can be 
computed from a uniform measure over blocked config- 
urations. The spatial correlation function obtained by 
a slow annealing close to p c is very similar to the one 
expected from the Edwards' measure. 

While these results and other presented below closely 
resembles the macroscopic behavior of many glassy sys- 
tems we will probe that this analogy can hardly be ex- 
tended to the organization of the phase space of the 
system and therefore suggests a non-universality of the 
glassy mechanism. 



III. RESULTS AND DISCUSSION 
A. The non-linear susceptibility 

It has been recently stressed that glassy dynamics has 
an inhomogeneous character. There one observes regions 
of correlated motion whose sizes depend of a control pa- 
rameter, and that survive for a characteristic time also 
parameter dependent jl^] . This character can be revealed 
studying the 4-point density correlation function. In lat- 
tice gas models this can be defined as 



where 



?(*) 



X (t)=N(( q \t))-(q(t)} 2 ) 



(1) 



(2) 



and the average is performed with respect to different 
samples. In p3[ , based on the theory of disordered sys- 
tems, it was proposed that in the deep supercooled phase, 
due to the inhomogeneous character of the dynamics, 
x(t) should display a maximum as a function of time. 
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This maximum appears more and more pronounced, and 
displaced to larger and larger times as the temperature 
is decreased (which is equivalent in the KA model to in- 
crease the density of particles in the system). 

This behavior reflects the geometrical structure of the 
phase space visited during the evolution of the system 
and the physics of time scale separation. Along its evolu- 
tion the system remains blocked in regions of phase space 
that have longer and longer time life, and with more and 
more correlated density fluctuations as the density is in- 
creased. These regions lye close to saddle points of the 
energy surface, which have lower and lower number of 
escaping directions, until at the dynamical transition the 
number of negative directions tend to zero. 

We have investigated the behavior of the function 
as a function of the density in the KA model, simulating 
systems with size L — 10 and L — 20, and averaging over 
100 samples. In figure ([!]) we show the behavior of the 
non linear susceptibility as a function of time for vari- 
ous values of the density, showing that qualitatively this 
function behaves in the same way found in simulations of 
liquids and in spin models (24) and theorized in |2^] . The 
maximum, as well as the position of the maximum seem 
to diverge at the dynamical critical density as power laws 
with exponents 4 and 2 respectively. 
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FIG. 1. The non linear susceptibility in the KA model as 
a function of time for different densities. From top to bottom 
p = 0.86, 0.85, 0.84, 0.83. The inset shows t max and \7yu 
close to the transition, the lines are power laws guide to the 
eyes with exponents 4 and 2 for t max and xTyn respectively 

It is interesting to note that in the low density phase, 
although the KA probes the configuration space ergodi- 
cally, the dynamics develops correlations which are sim- 
ilar to the ones seen in structural glasses. Even if the 
infinity time limit of the susceptibility is trivial, due to 
the fact that the system visit the space ergodically, the 
dynamics is highly correlated. 

In ref. |l8| ] it was emphasized that during the aging 
dynamics of the model, the measure concentrates on the 
blocked configurations. It remains to be explained why 



in the aging experiments the system blocks exactly at the 
value of the density where the fixed density experiments 
are not ergodic. 

B. The stability of the blocked states 

In order to understand the role of the blocked con- 
figurations as a function of the density, we studied the 
stability of the blocked states with respect to the move 
of a single particle, in a way analogous to study of the 
local stability in an energy landscape. 

To this end, we used the auxiliary model already intro- 
duced in jl8|. There, one associates to each configuration 
of the lattice gas an energy equal to the number of parti- 
cle which are free to move according to the KA dynam- 
ics. The thermodynamics of this auxiliary model can be 
studied easily numerically as a function of an auxiliary 
temperature T aux — l/(3 a ux- For example, the Edwards 
entropy as a function of the density can be obtained from 
the thermodynamic integration of: 

/>OC 

S(p) = H(p) + / d(3 aux (3 

JO 

where H{p) = —p log(/j) — (1 — p) log(l — p) is the entropy 
of the model at (3 aU x — 0. 

With the help of this auxiliary model we generated 
blocked configurations at different densities by lowering 
the auxiliary temperature slowly, in order to reach a typi- 
cal Edwards state at T aux = 0. We then moved a particle 
to an arbitrary site, and let the system evolve according 
to the KA dynamics. The perturbation to the blocked 
state causes some particle to become free to move, so we 
detected the number of particles that as a consequence of 
the perturbations move at least once. We observed two 
kind of events: 

• either the displacement of a particle has a little 
effect, and after a few particle move the system 
evolves again to a blocked state. 

• or the displacement of the particle completely un- 
blocks the system which starts to move ergodically, 
after a time which depends on the particle density. 

In figure || we display the estimated unblocking proba- 
bility, obtained generating 1000 blocked configurations 
using the simulated annealing as describe above and 
checking how many unblock, as a function of the density 
(white circles in figure 2). At low density, this proba- 
bility is different to zero in the low density phase which 
means that even if the system approaches by chance a 
blocked state, it will always depart from it. The value of 
the density at which the unblocking probability vanishes 
is compatible with the critical density p c = 0.88 of the 
dynamical instability. 
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We have also checked that configurations generated 
with a quenched of the auxiliary model form T aux = oo 
to T aux — give rise to the same unblocking probabil- 
ity, hinting that all blocked configurations have the same 
stability properties. 

0.9 n 1 1 1 1 1 1 1 n 
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FIG. 2. The unblocking probability as a function of the 
density, o Data from the block configurations obtained with 
the auxliary model, ■ Data from block configurations ob- 
tained taking a random configuration 



C. The geometrical structure of the blocked states 

As we have seen the blocked states are unstable for 
densities below threshold, while they are stable above 
threshold. One can ask if this stability property of the 
blocked states is in correspondence with their geometri- 
cal organization in space as it happens in mean field spin 
glasses p5fl . There the stationary points of the energy 
are saddles above threshold and minima below. It has 
been shown that the minima are disjoint below thresh- 
old: the maximal overlap between minima qo is smaller 
than one, and tends to one as the energy tends to the 
threshold energy from below, meaning that at that value 
of the energy one finds minima arbitrarily close from a 
given minimum. Above the threshold, there is no gap 
between saddles, and indeed one can show that the loga- 
rithm of the number of saddles at overlap q from a given 
saddle grows as q) « N(l—q). The threshold energy 
appears as a kind of percolation threshold for stationary 
point in the configuration space. It is worth noticing that 
while this organization is at the heart of the mechanism 
of dynamical transition, there is not sign in the Gibbs 
measure that reveals this structure. 

One might wonder if such an organization of sta- 
tionary points is also found in the KA model. In or- 
der to investigate this point we have studied the num- 
ber of blocked states at fixed overlap with a generic 
blocked states introducing in the auxiliary model a term 
N^—(q({ni},{n^}) — q) 2 , which for large e constrains 
the configurations to have the wanted overlap with the 



reference configuration. In the simulations the value 
e = 50 proved to be effective to get the wanted values 
of the overlap for all densities. We performed a simu- 
lated annealing of the model and we extracted the zero 
temperature entropy from thermodynamic integration of 
eq(l). 

In figure [| we show the entropy as a function of the 
overlap for different values of the density. It is apparent 
from the figure, that both above and below the transition, 
the curves behave qualitatively in the same way, with an 
entropy as a function of q which departs from zero reg- 
ularly at q = 1. This indicates that the structure of the 
blocked states does not change as we cross the transi- 
tion probing that the dynamic transition in this model is 
not linked to the geometrical organization of the blocked 
states. 
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FIG. 3. The entropy of the blocked states as a function of 
q for densities (from top to bottom) p — 0.8, 0.85, 0.9. 



IV. CONCLUSIONS 

In this paper we have investigated the mechanism 
which lead to vitrification in kinetic lattice gases, in com- 
parison with the one responsible of glassy phenomena in 
mean field models with random Hamiltonian. 

We have found that in the low density phase, despite 
the trivial character of density correlation implied by the 
equilibrium distribution, as it happens in supercooled liq- 
uids the dynamics has an heterogeneous character. The 
four point correlation function displays a maximum as a 
function of time that grows as a power of p — p c , and the 
dynamics remain correlated for times which also grow as 
a power of the same quantity. 

We have shown that the blocked states are unstable 
and play no role in the low density phase, while they 
are stable and act as attractor in the high density phase. 
Furthermore, we have investigated if this stability prop- 
erty corresponds to a different organization of the blocked 
states in phase space, with negative result. In both dy- 
namical phases (low and high density), one finds blocked 
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states at all distances from any blocked state. 

Then, we can conclude that despite the similarities in 
the dynamics, the glassy phenomena in kinetic models, at 
least at the level of the phase space organization, are very 
different from the ones found in systems with a complex 
energy landscape. 
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